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Abstract: A moisture detection of single rice grains using a slim and small open-ended coaxial 
probe is presented. The coaxial probe is suitable for the nondestructive measurement of 
moisture values in the rice grains ranging from from 9.5% to 26%. Empirical polynomial 
models are developed to predict the gravimetric moisture content of rice based on 
measured reflection coefficients using a vector network analyzer. The relationship between 
the reflection coefficient and relative permittivity were also created using a regression 
method and expressed in a polynomial model, whose model coefficients were obtained by 
fitting the data from Finite Element-based simulation. Besides, the designed single rice 
grain sample holder and experimental set-up were shown. The measurement of single rice 
grains in this study is more precise compared to the measurement in conventional bulk rice 
grains, as the random air gap present in the bulk rice grains is excluded. 

Keywords: small open-ended coaxial probe; single rice grain; measured reflection 
coefficient; relative permittivity; gravimetric moisture content; transmission line; 
microwave measurement techniques 



Sensors 2013, 13 



3653 



1. Introduction 

Recently, there has been an increased interest in microwave measurement techniques for the 
determination of moisture content, m.c. (or dielectric properties) of bulk rice grains in the microwave 
frequencies range [1-5]. This is due to the fact that the moisture percentage (wet basis m.c, ranging 
from 9% to 26%) in the rice grains plays an important role in the rice marketing and storage aspects. In 
marketing, the price of rice is dependent on the weight of the bulk rice, so accumulation of water in the 
rice grain will increase the price of rice. On the other hand, the moisture content in the grain mass 
determines the storage duration of the rice grains, in which m.c. < 13% indicates a storage duration of 
more than 60 days [6,7]. 

The m.c. of rice grains can be determined by either direct [8] or indirect methods [1—7]. Direct methods 
determine the m.cd by removing the moisture of the grains using oven drying methods (heating at 130 °C) 
or chemical reaction methods (extracting water using the reaction of iodine in sulfur dioxide). Both 
methods remove the moisture determine the true water content by the resulting weight loss. Indirect 
methods, in contrast, require the measurement of an electrical property of the grain using an 
instrument, a grain moisture meter. For very low frequency (DC) measurements, the desired electrical 
parameters are conductance, capacitance and resistance of the rice gain. The absorption power, resonant 
frequency, attenuation constant, reflection constant and transmission constant of the rice grain are the 
measured parameters of interest in microwave frequency measurements. Changes in electrical properties 
(or microwave properties) that can be directly correlated with a change in actual m.c. of the rice grain 
are obtained from the oven drying method (direct method). Recently, indirect methods have become 
more popular than direct methods due to their rapid tests and user friendly features. 

The prime considerations in the measurement of grain moisture using indirect methods are the size 
of the instrument sensor, which is in direct contact with the rice grain samples and the accuracy of the 
measurements [5]. Various microwave waveguide methods were proposed for the above purpose, but 
some of those methods require specify dimensions of the rice grain bulk to fit inside the given size of 
the waveguide [2-4]. However, the rice grain bulk is composed of a mixture of air and rice grains. The 
random distribution of rice grains and the air gap causes a low repeatability and low precision in the 
measurements. Among the mentioned methods, an open-ended waveguide method is the simplest and a 
nondestructive way to measure the m.c. of rice grains. The measurement using commercial open-ended 
waveguides is suitable for a specific rice grains size (width and length). This is because the waves 
scattered from the waveguide aperture would penetratd through small single rice grains as the small 
rice grains faed to entirely cover the aperture area of commercial probes. 

In this study, a millimeter size slim open-ended coaxial probe has been fabricated to measure the 
moisture in small grains. The coaxial probe was fabricated from a 0.86 mm outer diameter (OD) 
semi-rigid coaxial cable equipped with a male-type SMA plug connector. The coaxial cable was 
machined flat and polished to form an open surface end. Then, the coaxial probe was protected and 
covered by a customized stainless steel sheath with a flange, as shown in Figures 1 and 2. Here, the 
main focus is put on the aperture of the probe placed against a rice grain sample ranging from 9% m.c. 
to 26% m.c. Typically, the open-ended coaxial probe is calibrated by using open air, short terminator 
and liquid load. Directivity error, source match error, and frequency tracking are corrected by this 
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technique. Besides that, the rapid and simple calibration of the coaxial sensor without the use of short 
and load calibration kits were also proposed in this study (discussed in Section 3.2). 

2. Configuration and Dimensions Coaxial Sensor 

Figure 1 shows the side and the front sectional views of our milimeter size coaxial probe. The front 
sectional aperture of coaxial probe shows 2a = 0.20 mm diameter of inner conductor, 2b = 0.66 mm 
diameter of coaxial-filled Teflon and 0.86 mm diameter of the outer conductor. The coaxial-filled 
Teflon supports the coaxial line between the outer conductor and the inner conductor. Both inner and 
outer conductors guide the propagation wave in the coaxial line. In addition, an 11.3 mm total diameter 
steel flange is used to cover the total fringing field at the aperture probe. Figure 2 shows a picture and 
cross sectional structure of the coaxial probe. 

Figure 1. The side sectional view and front sectional view. 
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Figure 2. The Cross- sectional View. 
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3. Sensor Model 

3.1. Reflection Coefficient Model 

The reflection coefficient, r models of the millimeter coaxial probe are given as: 

r = &A7 3 +iM7 2 +2>X7 +2>X (1) 

n=0 n=0 n=0 n=0 

where the symbols S (in unit/ ), /? (in unit/ ), y (in unit/ ), and / are the complex coefficients for 
the polynomial Equation (1). The symbol/ represents the operation frequency. 

The complex parameters in Equation (1) was obtained by fitting the polynomial coefficients with 
calculated values obtained from the Finite Element Method using the COMCOL simulator over a 
broad range of permittivity values. The complex polynomial coefficients with seven decimals for 
Equation (1) are listed in Table 1. Equation (1) is valid for small coaxial probes, satisfying the relative 
permittivity, s r from 1 to 40 and the operation frequency from 0.4 GHz to 20 GHz. Comparison 
between the calculated and FEM simulated values for the reflection coefficient, T is shown in Table 2. 
If the FEM simulation results are used as the reference value, it is found that the percentage of relative 
error between both magnitudes of reflection coefficient will be less than 1%. 

Table 1. Complex Coefficient for Equation (1). 



Complex Coefficients in Equation (1) 



So 


-1.2489993 x 10" 42 
-j'2.8394016 x 10" 42 


Po 


4.068759 x 10" 32 
+ J6.0446093 x 10" 32 


7o 


-2.3422556 x 10" 22 
-yl. 0834847 x 10" 21 


Xo 


3.2710002 x 10" 13 
+ J4.3278983 x 10" 13 


d! 


2.8571222 x 10" 40 
+ j3 .768492 x 10" 40 


Pi 


-8.0370086 x 10" 30 
- /7.583469 x 10" 30 


?i 


4.4984192 x 10" 20 
+ jl. 5879853 x 10" 19 


Xi 


-6.0520681 x 10" u 
- ;5 .0636272 x 10" u 


8 2 


-2.3273679 x 10" 38 
-j 1.5437565 x 10~ 38 


P2 


5.7873753 x 10~ 28 
+ /2.8089148 x 10" 28 


?2 


-3.0910003 x 10" 18 
- j9.0522701 x 10" 18 


X2 


4.0103924 x 10" 9 
+ jl. 6441936 x 10" 9 


8 3 


7.6388153 x 10" 37 
+ j'6.4 15406 x 10" 38 


Ps 


-1.6900591 x 10" 26 
-y'8.6931643 x 10" 28 


73 


8.3505821 x 10" 17 
+ J2.6457801 x 10" 16 


X3 


-1.0502173 x 10" 7 
+ /7.2034269 x 10" 10 


84 


-6.5374999 x 10" 36 
+ j5 .4983527 x 10" 36 


P 4 


1.7398403 x 10" 25 
-j3. 1661421 x 10" 26 


74 


-7.5782519 x 10" 16 
- ;5.0073506 x 10" 15 


X4 


9.5824766 x 10" 7 
- ;2.7560669 x 10" 7 


8 S 


2.124032 x 10" 35 
- j 1. 2607055 x 10" 35 


Ps 


-4.0862657 x 10" 24 
+J1.4318571 x 10" 25 


7s 


3.2419268 x 10" 15 
+ J6.0109573 x 10" 14 


Xs 


-4.2081991 x 10" 6 
+ 338741 x 10" 6 


8 6 


-6.9741684 x 10" 35 
- j9.6327922 x 10" 36 


P6 


-2.122961 x 10" 25 
- J4.29076 x 10" 25 


76 


-8.983224 x 10" 15 
- J2.9066269 x 10" 12 


X6 


1.1443542 x 10" 5 
- j3 .7209467 x 10" 6 


8 7 


5.6295484 x 10" 35 
-j 1.6824543 x 10" 35 


P? 


-1.1659675 x 10" 25 
+ j'2.7664167 x 10" 25 


77 


7.1748551 x 10" 15 
- J7.8926933 x 10" 14 


X? 


9.9999087 x 10" 1 
+ ;'2.5383631 x 10" 6 
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Table 2. Calculated and Simulated Reflection Coefficient, T. 



f 

(GHz) 


Reflection Coefficient, r 
fe = l-;0) 


Relative 
Error (%) 


Reflection Coefficient, r 
fe.=40-y5) 


Relative 
Error (%) 


Equation (1) 


Simulations 


Equation (1) 


Simulations 


1 


0.9999955 
-j 0.002931223 


0.9999957 
-j 0.002926868 


0.00002 


0.9855305 
-j 0.09609176 


0.9824686 
-j 0.10045476 


0.2646 


10 


0.9995828 
-j 0.02934018 


0.9995701 
-j 0 .02929'/ '47 


0.00139 


0.5257350 
- j 0.7289402 


0.5155215 
- j 0.7370167 


0.0745 


18 


0.9986448 
- j 0.05295265 


0.9985924 
-j 0.05287694 


0.00563 


0.02905127 
- j 0.8579390 


0.02131934 
- j 0.85924107 


0.12505 



3.2. Calibration Model 



In this study; the simplest technique of de-embedding of coaxial probe is by extending the 
transmission phase in which the phase of reflection coefficient at measurement plane; AA' is extended 
towards the open end of coaxial probe; BB' using exponential term of exp (j2k c z). First; a full one-port 
calibration technique was implemented at the AA' plane using a commercial HP 85052D 3.5 mm 
calibration kit (open; short and load) which is only for network analyzer and cable error corrections. 
Secondly; under the assumption of quasi-TEM mode; the measured reflection coefficient; T^, of the 
sample at the plane AA' can be de-embedded to the end of the probe connector which coincide with the 
calibration plane BB' to give a reflection coefficient; T BB , by [9]: 



2jk c z 



(2) 



where z and k c = (2kJJc)^s c are the apparent physical length (in meter) and propagation constant of 
coaxial line, respectively. Symbols /, c and e c are the operation measurement frequency (in Hz), 
velocity of light in free space (299792458 ms" 1 ) and relative dielectric constant for the material filled 
in coaxial line (Teflon: e c = 2.05), respectively. In this work, only "open" standard involved in the 
calibration, the reflection coefficient measurements, Y Air for air at the plane AA' was taken, while the 
standard values for the air reflection coefficient, T Mr _ FEM at the plane BB' was simulated by Finite 
Element Method (COMSOL simulator). For the de-embedded process, the values of apparent physical 
length, z for coaxial line of the probe are required to determine. Once obtaining both values ( T^, and 



A_FEM 



), the apparent physical length, z can be calculated as: 



z = 



2k, 



In 



Air _FEM 



Air 



(3) 



= z+jz" 



Simultaneously, the attenuation constant, a in the coaxial line can be found from the optimized 
length, z as: 

(4) 



a = 





ft 

z 


I c J 


t 

z 
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Teflon and methanol liquid have been tested to verify the accuracy of the proposed calibration 
method. Figure 3 shows the comparison between the calibrated reflection coefficient measurement and 
the simulation results. In simulation, the relative permittivity, s r of Teflon was 2.05. While, the relative 
permittivity, s r , of methanol was computed by Cole-Cole model with parameters: s s = 33.7, e m = 4.45, 



x = 4.95 x 10 11 s and a = 0.036 [10]. The Cole-Cole model is: 



s = s„ + 



(5) 



\ + (jan) 



l-a 



In comparison, the calibrated measurement results using only the open standard have given a satisfy 
accuracy measurement up to 5 GHz. However, the noise inherent in both calibrated measurements are 
due to random errors contributed by the instruments or environmental measurement control. Moreover, 
the random errors are not taken into account in either calibration technique. The effect of the standing 
wave in the measurement of the open standard calibration becomes increasingly obvious when the 
operating frequency,/, is above 5 GHz. 

Figure 3. Variation in real and imaginary parts of reflection coefficient, (D\z (Tbb) and 3m 
(T BB )) at plane BB' with frequency, fat (25 + 1) °C. 
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The incident wave from the plane AA' is transmitted to the plane BB' by shifting phase of k c z, and is 
reflected back to input AA' with the same shifting phase. Thus, the aperture reflection coefficient, Y B b' 
at the plane BB' can be found by the phase delay of 2k c z with respect to the measured Taa' at the plane 
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AA' and the transmission line relationship is given as in Equation (2). However, the transmission line is 
imperfect and a fringing field occurs near the aperture probe. Hence, a phase shift between the forward 
wave and the reflected wave occurs, and produces the standing wave due to the superposition between 
the incident wave and the reflected wave inside the coaxial line [11]. The standing wave effect can be 
ignored if the operation quarter wave length, XIA in the coaxial line is large than the physical length, 
Z of coaxial line. For instance, 5 GHz of operation frequency will give XIA =15 mm which is smaller 
than the physical length, z ~ 22 mm, thus, the standing wave effect was significant when the operation 
frequency,/, was increased. 

3.3. Inverse Model 



For inverse solutions, the predicted values of the relative dielectric constant, s r , of a rice grain 
sample is obtained by minimizing the difference between the measured reflection coefficient, T B b' and 
Equation (1), Y by referring to the trial function, y/ : 

Y 3m(r) 3m(r„,)l (6) 



m 



Finding the zero routine was realized using the MATLAB fzero command. The initial approximate 
value in the numerical prediction was equal to 5 —j 0.001. The m terms in Equation (6) denote the 

3m(r BB ,) 

weighted parameter for ratio of j ^ . The weighted parameter, m is suggested due to the 

measurement with calibration using Equation (2) (only open standard: e,- = 1) does not consider the 
fringing effects for the medium loss samples. The predicted dielectric constant, s' r for Teflon and 

methanol liquid over frequencies 0.5 GHz to 12 GHz at room temperature (25 °C) are validated and 
compared with the FEM simulation, Cole-Cole models and the Agilent 85070E dielectric probe as 
shown in Figure 4. The diameter of outer radius for the Agilent 85070E probe (2b = 3 mm) is 
approximately 4.5 times larger than the studied coaxial probe. 

Figure 4. Variation in dielectric constant, s[ with frequency, / for Teflon and liquid 
methanol at (25 ± 1) °C. 
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4. Experimental 

The measurement reflection coefficient using millimeter coaxial probe that consists of the Agilent 
E5071C network analyzer in the frequency range between 0.5 GHz to 12 GHz was carried out at room 
temperature. The open end of the coaxial line was terminated by single rice grain sample. Normally, 
the length and width for various kinds of rice grains was in the range of 4.8-7.8 mm and 1.5-2.8 mm, 
respectively. Thus, the fringing field (sensing area ~ 2b) [9] from the coaxial probe aperture was 
sufficiently covered by the single rice sample. This can be applied based on the principle that the 
different signals are reflected from the terminal surface of the moist grain through the coaxial opening. 

Jati™ long grain white rice grown on the fertile soil of Kedah, the Rice Bowl of Malaysia, was 
used as the experimental sample. The rice grains were divided into different groups of 200 g per group. 
Each group of grains was sprayed with different estimated quantities of distilled water to achieve 
desired moisture levels. The bulk grains in each respective group were stirred and sealed in a container 
at 4 °C for 72 h to ensure a uniform water distribution within the bulk grains. The grains were 
conditioned to room temperature for 10 h prior to the measurements. Finally, 10 g of each group of 
bulk grain rice was dried in an air convection oven at 130 °C for 24 h [8]. The average moisture 
content, m.c. (in unit %) of each group of bulk rice grain was calculated on a wet basis as: 



(%) 



water 



xl00 



(7) 



m.c 



m 



waler 



+ m 



dry bulk grain 



where m 



■water 



and mdry bulk grain are mass of water and dry bulk grain, respectively. 



Figure 5. Experimental Set-up. 
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In this measurement, a specific holder was customized for the measurement of a selected single rice 
grain. The customized holder has a movable nylon platform which was mounted on a retort stand as 
shown in Figure 5. A single rice grain was randomly selected from each 10 g of bulk grain and placed 
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into a narrow and depth concave surface at the center of the platform. The steel flange of the coaxial probe 
was rigidly entered from the holder edge into a narrow space interval. The rice grain on the nylon platform 
was moved forward to the aperture coaxial probe by using a thumb screw. The concave circular surface 
on platform was used as a probe guide to ensure that the probe aperture is exactly touching the rice 
sample. The four springs were employed to ensure that the aperture probe contacts firmly with the 
interface single rice grain and to avoid the errors in the measurement of interface air gaps. 

5. Results and Discussion 

Figure 6 shows the reflection coefficient, Tbb; for 10 selected single rice grains from a bulk rice 
grain sample, where the bulk has a certain average water content, measured using the millimeter 
coaxial probe. As known, the m.c. for a bulk rice grain sample (which consists of thousands of single 
rice grains) is a statistical mean value due to the slightly different m.c S i ng ie of each single grain. Thus, 
the deviation of the actual m.c S mgU distribution in a single rice grain is higher as compared to the 
average m.c. of the entire bulk rice sample. This will lead to a scattered measurement of the reflection 
coefficient, Tbb; of the single grain in referring to the average m.c. of the bulk sample, as shown in 
Figure 6. However, Figure 6 does show a significant change of the measured reflection coefficient, 
T B B- of a single rice grain with the m.c. of bulk rice grain. The black solid line in Figure 6 is a 
regression fitting line from the average of measured reflection coefficient data (blue point-line). The 
average of the measured reflection coefficients was obtained from the 10 data points of the measured 
reflection coefficients for each bulk moisture content, m.c. The rice moisture calculations are based on 
a gravimetric method, thus differential effects of density between the bulk rice and the single rice grain 
in the measurement were ignored. In this study, the 2.44 GHz and the 5.81 GHz frequencies were 
chosen due to the fact those frequencies correspond to a free unlicensed band which is specifically for 
industrial, scientific and medical (ISM) measurement purposes. The 10.02 GHz band was also chosen, 
since the dielectric loss for the water is greater at around 10 GHz and thus provides a comparative 
approach for the rice measurements at higher frequencies. 

Figure 6. The relationship between reflection coefficient (magnitude, Tbb' and phase, 0) 
and the moisture content, m.c at 2.44 GHz, 5.81 GHz and 10.02 GHz, respectively. 
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Figure 6. Cont. 
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The polynomial regressions for the reflection coefficient, Tbb 1 (magnitude, \Tbb'\ and phase, 0) with 
respect to m.c. (in % units), are given in Equation (8) as listed in Table 3. Subsequently, Equation (8) was 
used by the coaxial probe to predict the m.c. (in % units), in a rice grain. It was found that, at higher 
frequencies, the point of measured reflection coefficient, T BB - shows a large variation of m.c. (in % units) of 
rice grains. 

The variations in relative dielectric constant, s' r , of rice grains with the percentage of m.c. at 2.44 GHz, 
5.81 GHz and 10.02 GHz, respectively, are plotted in Figure 7(a). The solid line of relative dielectric 
constant, s' r , in Figure 7(a) was the inverse of the reflection coefficient and refers to the trial function 
of Equation (6) with m = 1. The real part, 5He(Tsg') and imaginary part, JmiTss) in Equation (6) were 
calculated by using relationship of £He(rg 5 ')+/3m(rg B ')=|rgg'lexp(/0), where the values of |TWI and 0 
were obtained from Equation (8), while, the reference values of 9Kt(T) and 3m (T) in Equation (6) 
were computed using Equation (1). The measurement of dielectric constant points, e\ , of the single 
rice grain in Figure 7(a) were carried out using the studied coaxial probe with the Agilent 85070 
E.06.01.36 software in order to verify the data obtained from the dielectric inversion work. The 
tolerance, |A^'I of dielectric constant prediction between both techniques is shown in Figure 7(b), 
which gives the maximum deviation, \As' r | ~ 2 for all frequencies ranging 9.5% m.c. to 28% m.c . 
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Table 3. Polynomial functions for the reflection coefficient magnitude, \Tbb^, and phase, 0, 
(in rad units) respect to moisture content, m.c, of the bulk rice (in unit%) 

For/= 2.44 GHz, 

\T BB .\ = -4.918756 x 10" 5 m.c 2 + 3.694446 x 10" 4 m.c + 0.9948964, 
R 2 = 0.98420 (8a) 

0 = -3.043213 x 10~W 3 + 1.530463 x 10 3 m.c - 2.729901 x 10" 2 m.c + 0.1257916, 



R= 0.99338 



(8b) 



For /=5.81 GHz, 

|TWI = -1.616480 x 10~ 4 m.c 2 + 2.323330 x 10 3 m.c + 0.9810024, 
R 2 = 0.98386 (8c) 

0 = -6.330536 x 10~ 5 m.c 3 + 3.209160 x 10~W 2 - 5.643064 x 10" 2 m.c + 0.2498354, 



R 2 =0.98802 (8d) 



For / = 10.07 GHz, 

ITjB-1 = -2.135387 x 10~ 4 m.c 2 + 2.691858 x 10~ 3 m.c + 0.9783604, 
7? 2 = 0.98177 (8e) 

0 = -7.264877 x 10 5 m.c 3 + 3.722521 x 10~ 3 m.c 2 - 6.603589 x 10~ 2 m.c + 0.2683548, 

R 2 = 0.99068 (8f) 

• The symbol m.c refers to the percentage of bulk rice moisture content. 

Figure 7. (a) The variations in s' r of single rice grain with its bulk moisture content, m.c. 
at 2.44 GHz, 5.81 GHz and 10.02 GHz, respectively, (b) The absolute deviation, \&e' T I of 
dielectric constant prediction between the use of studied inversion technique and the 
Agilent 85070E software computation. 
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6. Conclusions 



The proposed coaxial probe has a small sensing area which covers the size of single rice grain and 
provides a nondestructive and real time moisture measurement for single rice grains. Moreover, the 
single grain measurement does not depend on the bulk density of the rice grains, thus the uncertainty 
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of bulk density in the rice measurement (due to different rates of broken rice in the bulk grain) can be 
ignored. In this study, moisture and dielectric models were created to suit the studied coaxial probe. 
The proposed simple de-embedding technique provides a rapid and low cost calibration procedure for 
the coaxial probe. However, the technique does not consider the systematic and random noises along 
the coaxial line, and thus it is suitable for a short coaxial probe, (z<A/4) due to the significant standing 
wave that is produced inside the long coaxial line. 
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